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The aim of this paper is to investigate the effect of a non-extensive form 
of statistical mechanics proposed by Tsallis on the formation of a quark- 
gluon plasma (QGP). We suggest to account for the effects of the dominant 
part of the long-range interactions among the constituents in the QGP by a 
change in the statistics of the system in this phase, and we study the rele- 
vance of this statistics for the phase transition. The results show that small 
deviations (~ 10%) from Boltzmann-Gibbs statistics in the QGP produce a 
noticeable change in the phase diagram, which can, in principle, be tested 
experimentally. 
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1. Introduction 



In the past few years the non-extensive form of statistical mechanics pro- 
posed by Tsallis [1] has found applications in astrophysical self-gravitating 
systems [2], solar neutrinos [3], high energy nuclear collisions [4], cosmic mi- 
crowave back ground radiation [5], and high temperature superconductivity 
[6, 7]. In these cases a small deviation of q (?» 10%) from Boltzmann-Gibbs 
(BG) statistics reduces the discrepancies between experimental data and the- 
oretical models. 

Substantial theoretical research has been carried out to study the phase 
transition between hadronic matter and the QGP [8-12]. When calculating 
the QGP signatures in relativistic nuclear collisions, the distribution func- 
tions of quarks and gluons are traditionally described by BG statistics. In 
this paper we investigate the effect of the non-extensive form of statisti- 
cal mechanics on the formation of a QGP. The crucial difference between 
the hadronic and QGP phase is the relative importance of short-range and 
long-range interactions among the constituents on either side of the antici- 
pated phase transition. The hadronic phase is characterized by a dominant 
short-range interaction among hadrons (which lends itself to BG statistics) 
while the QGP phase has a greatly reduced short-range interaction (due to 
"asymptotic freedom" ) and consequently a dominant long-range interaction. 
We suggest in this communication to account for the effects of the dom- 
inant part of this long-range interaction by a change in statistics for the 
constituents in the QGP phase. In addition we deem the non-dominant part 
of this long-range interaction negligible for the purpose of the phase diagram 
which we study here in detail. This latter view is supported by the empir- 
ical insensitivity of the phase diagram to details of the interaction among 
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constituents on either side of the phase transition. Therefore, we use the 
generalized statistics of Tsallis to describe the QGP phase while maintaining 
the usual BG statistics in the hadron phase (as we shall see this may also be 
regarded as choosing Tsallis statistics in the hadron phase with the Tsallis 
parameter q=l). 

Since hadron-hadron interactions are of short-range, BG statistics is suc- 
cessful in describing particle production ratios seen in relativistic heavy ion 
collisions below the phase transition [13-17]. Our motivation for the use of 
generalized statistics in the QGP phase lies in the necessity to include the 
long-range interactions on the QGP side. It has been demonstrated [18, 19] 
that the non-extensive statistics can be considered as the natural generaliza- 
tion of the extensive BG statistics in the presence of long-range interactions, 
long-range microscopic memory, or fractal space-time constraints. It was 
suggested in [4] that the extreme conditions of high density and temperature 
in ultrarelativistic heavy ion collisions can lead to memory effects and long- 
range colour interactions. Anticipating the formation of a QGP in heavy 
ion collisions, we use the generalized statistics to incorporate the effect of a 
long-range colour magnetic force in this phase. 

The generalized entropy proposed by Tsallis [1] takes the form: 

where k is a positive constant (from now on set equal to 1), w is the total 
number of microstates in the system, pi are the associated probabilities with 
J2iLiPi = 1) an d the Tsallis parameter (q) is a real number. It is straight- 
forward to verify that the usual BG logarithmic entropy, S — — J2T=iPi^ n Pi^ 
is recovered in the limit q — > 1. Only in this limit is the ensuing statistical 
mechanics extensive [1, 19, 20]. For general values of q, the measure S q is 
non-extensive. That is, the entropy of a composite system A © B consisting 
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of two subsystems A and B, which are statistically independent in the sense 
that pij 9 ^ = p\ A ^Pj B ^ , is not equal to the sum of the individual entropies as- 
sociated with each subsystem. Instead, the entropy of the composite system 
is given by Tsallis' q-additive relation [1], 

S q (A ®B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) (2) 

The quantity |1 — q\ can be regarded as a measure of the degree of non- 
extensivity exhibited by S q . 

The standard quantum mechanical distributions can be obtained from a 
maximum entropy principle based on the entropic measure [21, 22], 

S=-^[n i hin i ^(l±n i )Hl±ni)l ( 3 ) 

i 

where the upper and lower signs correspond to bosons and fermions, respec- 
tively, and fii denotes the number of particles in the i th energy level with 
energy e^. The extremization of the above measure under the constraints 
imposed by the total number of particles, 

E* = ^> ( 4 ) 

i 

and the total energy of the system, 

J2niei = E, (5) 

i 

leads to the standard quantum distributions, 

expf3(Ei - /i) =F 1 

where f3 — ^ and the upper and lower signs correspond to the Bose-Einstein 
and Fermi-Dirac distributions, respectively. 

To deal with non-extensive scenarios (characterized by q ^ 1), the ex- 
tended measure of entropy for fermions proposed in [6, 23] is: 

SPln,] = E{(^) + (7) 



which for q — > 1 reduces to the entropic functional (3) (with lower signs). 
The constraints 

E"' = N (8) 

i 

and 

= £ (9) 

i 

lead to 

Wi = j (10) 

[1 + ^-1)^- ^5=1+1 

In the limit q — > 1 one recovers the usual Fermi-Dirac distribution (6) (with 
lower sign). 
Similarly, 

n l = 1 (11) 

[1 + ^-1)^-^-1 

for bosons. We now turn to the description of the system in the QGP phase. 



2. The QGP phase 

In this approach, the quarks and gluons are treated as forming an ideal 
gas, apart from the non-perturbative corrections to the pressure and energy 
density resulting from the bag model [8]. We initially (a) describe the system 
by BG statistics in order to define our notation and general numerical pro- 
cedure. In (b) we lay out the differences due to the non-extensive statistics. 
We emphasize that only (b) can incorporate the anticipated long-range forces 
in the QGP phase. 
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(a) Boltzmann-Gibbs (BG) Statistics 



According to the BG statistics the energy density, pressure and baryon 
number density for a QGP consisting of massless u and d quarks and anti- 
quarks, each with degeneracy factor cIq = 12, and gluons, with degeneracy 
do = 16, at a temperature T and baryon chemical potential \i are given by 

do r r°° ,, r _ _ , d ( 



U QGP = ^{ / dk k 3 [n Q + n Q ] + ^ J dk k 3 n G + B (12) 

d roo I _l 

Pqgp = 2^2" J o dk k 2 {ln[l + exp -(/i Q - k)} + ln[l + exp —(/iq + k)} 

/ dk k 2 ln[l - exp(— )] - B (13) 
,/o T 



and 



where 



2vr 2 



n QGP 



roo 

{ dkk 2 [n Q -n Q ] (14) 

J 



6tt 2 

™Q(<5) = — m.t .. n i 1 ' ( 15 ) 



n G = ZTTTT^ 7 ( 16 ) 



exp^(A;=F^Q) + 1' 
1 

exp(|) - 1 
and fi Q = |. 

Integration by parts of equation (13) yields 

Pqgp = ^(uqgp ~ 45) (17) 

Integration of equations (12)-(14) yields 

_ * 2 (f1 | 7 . ^ T 4 | d Q fi 2 T 2 d Qf i 4 
uqgp - ^(d G + -d Q )T + + ^- 2 +B (18) 

_ * 2 ( j | 7 , )T 4 . dQfi 2 T 2 d Q /x A 
PQGP - ^(d G + ~ A d Q )T + + - 5 (19) 



and 



n QCF = d Q [^ + -Lj (20) 
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where B is the bag constant which is taken here as (210 MeV) 4 [9] with an 
uncertainty of ~ 15%. 



(b) Generalized Statistics 

If we use the generalized statistics to describe the entropic measure of the 
whole system, the distribution function can not , in general, be reduced to a 
finite, closed, analytical expression [23-28] . For this reason we use generalized 
statistics to describe the entropies of the individual particles, rather than of 
the system wholef]. 

The distribution function is then given by 

n t = 1 , (21) 

[l + ^-l^-^FT^l 

where the upper and lower signs correspond to bosons and fermions, respec- 
tively [23]. In this case 

1 

[l + ±(g-l)(fc=F^)]^ + l 
and 

n G = ■ 1 (23) 

[l + i(g-l)fc]i=r-l 

In the limit q — > 1 one recovers equations (15) and (16). 
The expression for the pressure is given by 

d Q T r,,,,,/;v ~> , /r-i. 



n Q(Q) = i (22) 



PQGP=^ dkk^^— + ^—)- 

27r jo q — 1 q — 1 

dGl f dkk 2 (^ -)-B (24) 



2tt 2 Jo q 



1 For a more detailed account of this important point see ref. [23]. 
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where 

f Q = l + [l + ^(q-l)(k-» Q )}^ (25) 
f Q = l + [l + ^(q-l)(k + iiQ)]^ (26) 

and 

fa = l-[l + ^(q-l)k]^ (27) 
which in the limit q — > 1 reduces to (13). 

Since the integrals in (12)-(14) are not integrable analytically one has to 
calculate these integrals numerically. For q > 1, the quantity [1 + ^(g — 
l)(fc — /xq)] becomes negative if /xq > k. To avoid this problem we use [29], 

f Q = l + [l + ^(q-l)(k-fi Q )}^, k>^Q (28) 

and 

fq = 1 + [1 + ^(1 -q)(k-fJL Q )]^, k<^ Q (29) 

In the limit q — > 1 one recovers, of course, the appropriate Fermi-Dirac 
distribution in both cases. We now turn to the hadronic phase of the system 
which is more readily accessible to experiment. 

3. The hadron phase 

The hadron phase is taken to contain only interacting nucleons and antin- 
ucleons and an ideal gas of massless pions motivated by the findings in [30] . 
The interactions between nucleons can be treated either by means of an 
excluded volume approximation or by a mean field approximation. The ex- 
cluded volume approximation to hadron interactions is either thermodynam- 
ically inconsistent, or leads to computational difficulties in realistic calcula- 
tions if the consistency problem is addressed. To treat the interactions we 
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use the relativistic mean field or Hartree approximation proposed by Walecka 
[31] which is thermodynamically self-consistent. In this model the interaction 
between nucleons is described by the scalar-isoscalar a and vector-isoscalar 
uj^ mesonic fields with baryon-meson interaction terms in the Lagrangian: 
g^ipa and g^tp^ipu^. For nuclear matter in thermodynamical equilibrium 
these mesonic fields (a and u^) are considered to be constant classical quanti- 
ties. The scalar field a describes the attraction between nucleons and lowers 
the nucleon (antinucleon) mass M to M* = M — g a (a) . The nucleon-nucleon 
repulsion is described by the vector field u 11 which changes the nucleon (antin- 
ucleon) energy by (±U(n)) Q. 

The thermodynamically self-consistent equations of state (EOS) for nu- 
clear matter are [32, 33]: 

p( t , A*) = -7T / t \z r== ^( njv + + nU 
3 J (2?r) yk 2 + M* 

rn 

/ dn'UM) + P(M*) (30) 
J o 



n N (N) = [exp(- ) + 1] (31) 

5P , dP{M*) r d 3 k M* . n 

)t,„ = -In I /0 _ v , ; = (n N + %) = (32) 



<5M* dM* MV 7 (2tt) 3 /p ~ M *2 



V)=lN J J^^( n N ~ u n) (33) 



and 



//73 L I „ n 
— ^k2 + M* 2 (n N + n R )+j o dn'U(n')-P(M*), (34) 

where p, n, and u are the pressure, baryon number density, and energy density 

respectively, Un(n) is t ne distribution function of nucleons (antinucleons), fj, 
2 The odd G-parity of the w-meson is responsible for the attractive w-exchange in NN 
scattering as compared to the repulsive w-exchange in NN(NN) scattering. 
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is the baryon chemical potential, and 7jv is the spin-isospin degeneracy of the 
nucleon which is 4 for symmetric nuclear matter. Equation (32) describes 
the dependence of the effective nuclear mass M* on T and /i which is defined 
by extremizing the thermodynamical potential (maximizing the pressure). 
If we choose [31], 

P(M>-1(M-Mf, U{n) = Cln (35) 

where C a = g a (^), and C w = g u {^)- The parameter set (M, g a , g u , m a , 
mj) consists of phenomenological constants determined by experiment. We 
take here M= 0.940 GeV, M*= 0.543M, g a « 11 GeV" 1 , g u « 14 GeV -1 , 
m a = 0.520 GeV, m u = 0.783 GeV, which reproduces data (see ref. [33, 34] 
for details) in the hadronic phase. 

The energy density, pressure, and baryon number density for the hadron 
gas taken to contain nucleons, antinucleons and an ideal gas of massless pions 
are: 

u H (T, //) = £L J°° dk k 2 yjk* + M* 2 (n N + n R ) + ^C> 2 + 

Z7T JO 2 

1 ; (M-M*) 2 + — tt 2 T 4 (36) 



2cr io 



lrn , 7at f°° dkk^ , , 1 2 

p H (T,fi) = — ( nN + nx) + -C u n 

671 Jo ^Jk 2 + M* 2 



' (M - M*f + — 7r 2 T 4 (37) 



2C^ ' 30 

and 

/■y^ roc 

n H {T,fi) = — dkk(n N -n N ), (38) 
Itx a Jo 

where fi N ^ as in (31) with U(n) as in (35). 

Finally, we address the phase transition from the hadronic to the QGP 
phase, within our model. 
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4. Phase transition 



Assuming a first order phase transition between hadronic matter and 
QGP one matches an EOS for the hadronic system and the QGP via Gibbs 
conditions for phase equilibrium: 

Ph = Pqgp, T h = Tq GP , fi H = /j,q G p (39) 

With these conditions the pertinent regions of temperature T and baryon 
chemical potential fi are shown in fig. 1 for q= 1, 1.1 (0.9) and 1.25 (0.75). 
The critical temperature at \x = for q= 1,1.1 (0.9) and 1.25 (0.75) are found 
to be 148 MeV, 122 MeV and 79 MeV, respectively. As the non-extensive 
parameter deviates from q= 1 to 1.25 (0.75), the critical temperature becomes 
almost independent of the baryon chemical potential which is associated 
with the number of particles. The variation of the bag constant B between 
(180 MeV) 4 and (250 MeV) 4 does not alter our findings significantly. One 
still observes a flattening of the T(fj) curves in fig. 1 as |1 — q\ increases. 
The only effect is a shift of the value of the maximal ji and T(0) in fig. 
I. There is no significant change of the slope of T(/x) over the depicted 
range. The agreement between BG and the generalized statistics as the 
critical temperature approaches zero is evident from equations (13) and (24). 
Fig. 2 shows the dependence of the critical temperature on the non-extensive 
parameter at fj, — 0, 250 MeV and 1000 MeV. The dependence is almost linear 
for small values of q with a slope |^-| ~ 240 MeV. This can be interpreted as 
a new type of universality condition [6] which suggests that the formation of 
a QGP occurs (almost) independent of the total number of baryons in heavy 
ion collisions. 
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5. Conclusion 



We have studied the phase transition from a system in the hadronic 
phase to the QGP phase. On the hadronic side the detailed form of the 
(short-range) interactions among hadrons turns out to be unimportant for 
the phase transition. On the QGP side the short-range interactions have 
died out (due to "asymptotic freedom") and have made room for the only 
remaining long-range interactions among the constituents. We take account 
of the dominant part of this interaction by a change in the statistics of the 
system in the QGP phase. We present here testable consequences of us- 
ing the non-extensive form of statistical mechanics proposed by Tsallis in 
the QGP. The resulting insensitivity of the critical temperature to the total 
number of baryons presents a clear experimental signature for the existence 
of non-extensive statistics for the constituents of the QGP. 
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Figure 1: Phase transition curves between the hadronic matter and QGP for 
q=l (solid line), q=l.l (0.9) (dotted line) and q=1.25 (0.75) (dashed line). 
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Figure 2: The dependence of T on |1 — q\ at /i— (solid line), /i— 250 MeV 
(dotted line) and fjt= 1000 MeV (dashed line). 
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